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ABSTRACT 

In this paper we explore collisional ionization and time-evolving photoion¬ 
ization in the. X-ray discovered, ionized absorbers in Seyfert galaxies. These 
absorbers show temporal changes inconsistent with simple equilibrium models. 

We develop a simple code to follow the temporal evolution of non-equilibrium 
photoionized gas. As a result several effects appear that are easily observable; 
and which, in fact, may explain otherwise paradoxical behavior. 

Specihcally we hnd that: 

1. In many important astrophysical conditions (OVII, OVIII dominant; and 
high lO^^'^cm”^) column density) pure collisional and photoionization 
equilibria can be distinguished with moderate spectral resolution observa¬ 
tions, due to a strong absorption structure between 1 and 3 keV. This feature 
is due mainly to iron L XVII-XIX and Neon K IX-X absorption, which is 
much stronger in collisional models. This absorption structure may be mis¬ 
interpreted as a flattening of the intrinsic emission spectrum above ~ 1 keV 
in low resolution data. 

2. In time-evolving non-equilibrium photoionization models the response of the 
ionization state of the gas to sudden changes of the ionizing continuum is 
smoothed and delayed at low gas densities (usually up to 10® cm“®), even 
when the luminosity increases. The recombination time can be much longer 
(up to orders of magnitude) than the photoionization timescale. Hence a 
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photoionized absorber subject to frequent, quick, and consistent changes of 
ionizing luminosity is likely to be overionized with respect to the equilibrium 
ionization state. 

3. If the changes of the ionizing luminosity are not instantaneous, and the 
electron density is low enough (the limit depends on the average ionization 
state of the gas, but is usually ~ 10^ cm“^ to ~ 10® cm“®), the ionization 
state of the gas can continue to increase while the source luminosity de¬ 
creases, so a maximum in the ionization state of a given element may occur 
during a minimum of the ionizing intensity (the opposite of the prediction 
of equilibrium models). 

4. Different ions of different elements reach their equilibrium conhguration on 
different time-scales, so models in which all ions of all elements are in pho¬ 
toionization equilibrium so often fail to describe AGN spectral evolution. 

These properties are similar to those seen in several ionized absorbers in 
AGN, properties which had hitherto been puzzling. We applied these models to 
a high S/N ROSAT PSPG observation of the Seyfert 1 galaxy NGG 4051. The 
compressed dynamical range of variation of the ionization parameter U and the 
ionization delays seen in the ROSAT observations of NGG 4051 may be simply 
explained by non-equilibrium photoionization model, giving well constrained pa¬ 
rameters: He = (l.O^oJ) X 10® cm“®, and R= (O.Td+o;®))) x 10^® cm (~ 3 light 
days). 
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1. INTRODUCTION 


The detection and the study of ionized absorbers is more difficult than that of cold 
neutral absorbers, but can yield much more detailed information about the nature of the 
absorbers and the state and geometry of the nuclear regions of AGN. If photoionization 
applies then the electron density of the gas and its distance from the ionizing source can be 
estimated. 

Absorption features from ionized gas, are common in the X-ray spectra of Seyfert Galax¬ 
ies and some quasars Q Deep oxygen VII and VIII absorption edges at 0.74 keV and 0.87 
keV (rest) have been detected by the ROSAT PSPC and the ASCA SIS in quite a large 
number of Seyfert 1 galaxies (Reynolds, 1997). 

We do not know yet what is the origin of the gas ionization. Models to date have 
assumed the simplest equilibrium photoionization case (e.g. Fiore et ah, 1993; Guainazzi 
et ah, 1996; Reynolds et ah, 1995). In this case, if the gas is conhned in a single cloud 
of constant density, and if the recombination time is smaller than the typical variability 
timescale, then the ionization parameter, and hence the ionization state of the gas, should 
follow closely the intensity of the ionizing continuum. This is not always observed. In 
two ASCA observations of MGC-6-30-15 (Fabian et ah 1994, Reynolds et al. 1995) the 
best fit ionization parameter is higher when the ionizing flux is lower, in contrast with the 
expectations of the simplest equilibrium photoionization model; the ASCA observations of 
MR2251-178 show a roughly constant ionization parameter despite large variation in the 
2-10 keV flux (Reynolds & Fabian, 1995); finally, in a ROSAT observation of NGC 4051 
(McHardy et ah, 1995), the ionization parameter does not linearly track the luminosity, but 
shows changes that are smoothed and delayed with respect to the luminosity changes. It 
seems clear that, at least in the above three cases (which are also the best studied), the 
simplest photoionization equilibrium model is inadequate. We clearly need more complete 
and consistent models to interpret the available data. The gas could well be distributed 
in an irregular region with varying density. Different ionization states would then apply to 
different region of this gas. Other authors have adopted such multi-zone models (Otani et 
al. 1996, Kriss et al. 1996). In those models the authors assume that absorption features 


^ NGC 4051: Me Hardy et al. 1995; Guainazzi et al., 1996; MGC-6-30-15: Nandra & Pounds, 1992, 
Reynolds et al. 1995, Orr et al., 1997; NGC 3783: Turner et al. 1993, George et al., 1996, George et al., 
1998; NGC 985, Brandt et al. 1994, Nicastro et al., 1997; NGC 5548: Done et al. 1995, Mathur et al. 1995; 
NGC 3227, Ptak et al. 1994; NGC 3516, Kriss et al. 1996; IC 4329A, Cappi et al. 1996; PG 1114-^445: 
Laor et al. (1997), George et al., 1997, Mathur et al., 1998; MR 2251 — 178, Halpern (1984), Pan, Stewart 
and Pounds (1990); 3C 351, Fiore et al. (1993), Nicastro et al., 1998; see also Reynolds (1997) 
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from different ions of the same element are imprinted on the spectrum by the transmission 
of the ionizing radiation through multiple distinct clouds of gas with different geometrical 
conhgurations, ionization states and densities. This is of course a possibility, but to us it 
seems rather ‘ad hoc’. 

In this paper instead we discuss the additional physics of collisional ionization, and of 
time-evolving photoionization in a single zone model. Time-evolving photoionization models 
have been applied in the past to the extended narrow line emitting region of radio galax¬ 
ies, to predict the evolution of the line ratios in low-density (ue = 1 cm“^), low ionization 
(LogU r'U -3), multi-zone photoionized gas after the switch-off of the ionizing source (Binette 
et ah, 1987). We present here an iterative solution of the time-dependent ionization balance 
equations (accounting for all ions of the most abundant elements in gas with solar composi¬ 
tion), and apply our models to a drastically different physical and geometrical conhguration: 
a relatively high-density {ue = 10® — 10® cm“®), medium-high gas ionization (LogU ~ -1, 
1), geometrically thin {DeltaR/R ~ 10“^), single zone absorber along the line of sight, 
undergoing rapid and persistent variations of the ionizing flux. 

As an example, we apply our models to the ROSAT data of the Seyfert 1 galaxy 
NGC 4051. A detailed analyses of the time behaviour of ionized gas in Seyfert 1 galax¬ 
ies, using the ASCA-SIS relatively high resolution data, is deferred to a future pubblication. 


2. Ionization Models 


The innermost regions of an AGN are likely to be ‘active’, in the sense that the gas 
there conhned is expected to be involved in signihcant bulk motion, as in strong outflows 
(Arav et ah, 1995) or inflows. In particular the evidence for ionized outflows from the cen¬ 
tral regions of AGN is strong (Mathur et al. 1994, 1995, 1997, 1998). Hence sources of 
mechanical heating of the gas (for example adiabatic compression by shock waves) may well 
be at work in the high density clouds, making collisional ionization the dominant ionization 
mechanism. In the low density clouds the gas could be far from photoionization equilib¬ 
rium, because the recombination timescales have become greater than the X-ray continuum 
variability timescales. Hereinafter, by photoionization time (tph) and recombination time 
{tree) we mean the time necessary for the gas to reach photoionization equilibrium with 
the ionizing luminosity during increasing and decreasing phases respectively. We discuss in 
turn the possible importance of collisional ionization and time-evolving photoionization, in 
determining the final transmitted spectrum that we observe in AGN. 



2.1. Collisional Ionization Models 


To study collisional ionization we constructed a series of pure collisional ionization 
models and compared them with the equilibrium photoionization models. The models were 
created using CLOUDY (version 90.01, Ferland, 1996), fixing the ionization parameter (i.e. 
the dimensionless ratio between the number of Hydrogen ionizing photons and the electron 
density of the gas) U to 10“^ and calculating the spectra transmitted through clouds with a 
total column density, Nh = 10^^'® cm“^, constant density and temperatnre (the calculation 
of the physical condition of clouds in coronal or collisional equilibrium is one of the options of 
CLOUDY; Ferland, 1996). The distribntion of fractional abundances for the most important 
ions is essentially independent of the electron density valne: for 10 < Ue < 10^^ cm“^; the 
variation is smaller than 13% for ions with fractional abundances greater than 0.01. 

For temperatures Tg in the range 10^ — 10^'^ the main edges imprinted on the spectra 
are ronghly the same as those in the spectra from photoionized gas with —1 < log{U) < 1.5. 
However, the distribution of the different ionic species is very different in the two cases, and 
so the relative optical depth of the edges differ markedly. High quality X-ray spectra, where 
more than one edge is visible, could then discriminate between the two cases. 

The different ionic abundances in the collisional and photoionized cases can be seen in 
Figs. la,b and 2a,b. Here we plot the fractional abnndances of NeVHI-NeXI and OVI-OIX 
(Fig. 1) and FeXIV-FeXXV (Fig. 2), computed in the case of photoionization and collisional 
ionization respectively as a fnnction of U and Tg. As Fiore et al (1993) and Mathur et al. 
(1994) have shown, the photo-ionization state of the gas depends strongly on the spectral 
energy distribntion (SED) of the ionizing continnum, from radio to hard X-rays. To simplify 
the comparison between photoionization models and the data, we therefore used an ionizing 
continuum similar to the SED of NGC 4051 (Done et al 1990). 

Firstly we note that the fractional abundances in photoionized gas are more smoothly 
distribnted than those of collisionally ionized gas. Let us consider the regime in which 
the OVH and OVHI fractions are higher than si 0.2, so that both OVH and OVHI edges 
are present in the emerging spectrnm, as sometimes fonnd in Seyfert galaxies spectra (e.g. 
Otani et al. 1996, Gnainazzi et al. 1996). In the photoionization case the range of U where 
this occurs is 1.2-3.6, a factor 3. This is twice as wide as the corresponding range of Tg: 
1.5 X 10® — 2.4 X 10® K (a factor of 1.6: dotted lines in panels a and b of Fig. 1). Moreover, 
NeIX is abnndant and Garbon is almost fully stripped in the collisional case, bnt not in the 
photoionized case. 

The difference in the resnlting spectra is shown in Fig. 3a. Here we show the ratio 
between two power law spectra emerging respectively from clouds of collisionallly ionized 
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(Te = 1.8 X 10® K) gas, and photoionized (U=2) gas with solar abundances (the choices of 
Te and U were made by looking for similar OVII-OVIII abundances: see Fig. la). Below 
1 keV there is a large feature due to the deep CVI K absorption edge in the photoionization 
case. The most important feature above 1 keV is the NeIX absorption edge in the collisional 
case, but this feature is smoothed by the presence, in both cases, of deep OVII and OVIII 
absorption edges. Reynolds (1997), htted both a simple OVII and OVIII K-edges model 
and a physical photoionization model to ASCA SIS spectra of a sample of Seyfert 1 known 
to host a warm absorber. Comparing the results he found that in some cases, the best 
fit 2-edges model continue to show residuals at the energy of CVI K edge, while residuals 
to physical photoionization model do not. This may point to photoionization as ionizing 
mechanism. Spectra with higher resolution from 0.4 keV to 3 keV should therefore allow one 
to diagnose absorption by collisionally ionized or photoionized gas that has significant OVII 
and/or OVIII abundances (e.g. Otani et al. 1996, Guainazzi et ah 1996). 

The other important feature of collisional models is the inertia of the heavy elements 
(from Ne to Fe) against becoming highly ionized, even when the oxygen is almost fully 
stripped. In Fig. la and lb we have marked the range of U for which OIX is by far the 
dominant ion (no/x~0.8) and OVIII is the only oxygen edge in the emergent spectrum 
(r^O.8, for an of 10^^'® cm“^ and solar abundances). These same intervals of U and Te 

are also shown in Fig 2a and 2b. In the photoionization case Neon is almost fully stripped 
(NeXI), and the dominant ions of the iron are FeXVIII to FeXXIII. Instead in the collisional 
case there is a range of temperature for which Neon and Iron are much less ionized, with 
the dominant ions being NeIX-NeX, and FeXVI to FeXX. In particular we note the large 
abundance of FeXVII (FeXVII and FeXXV are respectively Ne-like and He-like, and so very 
stable). The different ionization level of Neon and Iron implies a very different emergent 
spectrum between 1 and 3 keV. 

This is illustrated in Fig. 3b which shows the ratio between spectra transmitted from 
collisionally ionized gas with Te = 3.9 x 10® K, and photoionized gas with U=10 (both with 
Nh = 10^^ ® cm“^). With these values of Te and U, noix ~ 0.8 in both collisional ionization 
and photoionization models (see Fig. la). The ratio does not show any signihcant feature at 
E< 1 keV (implying similar ionization states of Carbon, Nitrogen and Oxygen), but between 
1 and 3 keV, the spectrum from collisionally ionized gas shows a large and complicated ab¬ 
sorption structure due mainly to iron L XVII-XIX and Neon K IX-X absorption (Fig. 3b). 
This absorption structure may be mis-interpreted in moderate quality spectra as a flatten¬ 
ing above ~ 1 keV. Again, higher spectral resolution data with good S/N can distinguish 
between photoionization or collisional equilibrium. 

That the NeIX and FeXVII edges are of similar depth is due to a coincidence of cross- 
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sections and abundances. The photoelectric cross section of the iron L shell atoms ranges 
from ~ 2 X cm^ for the Lp levels of FeXVI-FeXIX {E - (1.17-1.47) keV), to ~ 3 x IQ-^o 
cm^ for the Lp level of FeXXII (E=1.78 keV) (Kallman & Krolik, 1995). 0 Furthermore 
in the collisional case the relative abundance of FeXVII is in the range 0.2-0.7, compared 
to ;;0.25 in the analogous photoionization case. The relative abundance of the NeX, in the 
collisional case, is greater than 0.25, but NeX ^0.4 in the analogous photoionization case, 
while 0.051; (NelX)coZZ 10.40 and (NelXjphot 10.05 (Fig. 1 a,b). The K-edge energy of NeX 
is 1.36 keV the same of that of the FeXVIII Lp-edge. While the photoelectric cross-section 
of Neon is about a factor 4 lower than Iron, the solar abundance of the Neon is about a 
factor 4 greater than that of Iron (Grevesse & Anders 1989). Furthermore the cross-section 
of NeIX is about a factor 3 higher than NeX. As a result, spectra from a cloud of collisional 
ionized gas at Tg ~ 4 — 6 x 10® K, will show a very deep (T^e^pxy// ~ 1.0 x edge 

at the FeXVII Lp-edge energy (1.26 keV), a deep {rp^Lpxviii+NeKx ~ ^ ^H 2 i) edge at 

1.36 keV due to the comparable contributions of the NeX K-edge and the FeXVIII Lp-edge, 
and a similarly deep absorption edge at the NeIX K-edge energy (NeIX) = 1.196 keV, 
Tne^ix ~ 0-6 X Vhjs)- In the corresponding photoionization case, the relative abundances 
of both NeIX-X and FeXVII-XVIII are too low (Fig. la,b, 2a,b) to imprint similar features 
on the spectrum. 

When collisions are the dominant ionization process, ([/;;0.01) we can estimate the 
minimum distance, R, between the ionized gas cloud and the central source. For an ionizing 
luminosity Ljon, we hud: i? > 1.6 x 10 ^®(neio)“®'®(Q 52 )®'® (with rZeio being the electron density 
in units of 10^® cm“®, and Q 52 the rate of photons ionizing hydrogen in units of 10®^ ph s“^: 
Q ~ I^Lyd-E{Lion{E)/E)). For typical AGN ionizing continuum shapes and luminosities 
{Lion = 10"^^, 10"^® erg s“^), and assuming = 1, we hnd R> 6.2 x 10^® cm, 2 x 10^® cm, 
respectively, similar to the size of the BLR for such AGN (Peterson et ah, 1993). 


2.2. Mixed Collisional and Photoionization Models 


The above discussion concerns gas in pure collisional equilibrium. If the gas density is 
low enough and/or if the gas is close enough to the central X-ray source then photoionization 
can be important too. We examined a number of models with a varying mixture of collisional 


^The photoelectric cross section of the Ls levels of the iron ions, is about one order of magnitude lower 
than the correponding Lp cross sections, and the Lp levels of the ions FeXXIII-FeXXVI are not populated 
in ordinary conditions. 



ionization and photoionization. For T^. = 2.5 x 10® we find that for values of U ^0.1 the main 
features of the fractional abundance distribution are still those of pure collisional ionized gas. 
Increasing U to 0.3 causes a sharp change in the ionization structure. For 0.3 < U ^100 
(the exact value of the upper limit depending on the equilibrium temperature determined 
by collisional mechanisms) the ionization structure is determined by both processes, and the 
transmitted spectra resemble those from pure photoionized gas of much higher ionization 
parameter. This is shown in Fig. 4 where we plot in the lower panel three spectra from 
purely photoionized gas (column density logiVj7=22.5 and log(U)=0.7, 1.1 and 1.4), and in 
the upper panel 5 spectra from gas in which both mechanisms are at work (logiVj:^=22.5, 
Te = 2.5 X 10® and log(U)=-2, -1, 0.3, 0, 0.4). The pure photoionization spectra with 
log(U)=0.7, 1 (lower panel) are very similar to the mixed collisional and photoionization 
spectra with a factor 4-10 times lower U (log(U)=-0.3, 0.4; upper panel). Fitting the mixed 
spectra with pure photoionation models would give a much more compressed range of U 
than the real one, as found by McHardy et ah in NGC 4051. The above spectra in the two 
panels of Fig. 5 are practically undistinguishable even at high resolution, because the shape 
of the fractional abundance distributions when both processes are working, closely resemble 
those of a purely photoionized gas (see upper panels of Fig. 1 and 2). Fortunately the two 
models predict quite different delay properties (§2.3). 


2.3. Time-evolving Photoionization Models 


If photoionization is the dominant process (the gas being purely photoionized, or the 
density enough low to give U70.1) the main features of the fractional abundance distribution 
are those of photoionized gas. However, if the X-ray source is variable, photoionization 
equilibrium will apply only if the density is high enough to make the ion recombination 
timescales shorter than the variability timescales. There are regimes where photoionization 
equilibrium does not apply at all. Here we investigate the low density case in some detail, 
beginning with a discussion of the relevant physics. 


2.3.1. Equilibrium photoionization 


Let us suppose that a geometrically thin single cloud of optically thin, gas is illuminated 
by an intense flux of ionizing radiation emitted by a variable source located at a distance R 
from the cloud. 
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We can calculate the equilibrium distibution of the ionic species in the cloud, by setting 
the photoionization rate equal to the radiative recombination rate | (see, e.g., Netzer, 1990). 


('X') f 

V nxi J arec{X\Te)ne 

and adding the condition for charge conservation = 1- Here arec{.^\Te) is the 

radiative recombination coefficient (in cm^ s“^), which includes recombination to all levels, 
and Fxi is the photoionization rate of the ion X*, which, for optically thin gas clouds, can 
be written (see, eg., Netzer, 1990): 


( 2 ) 


Fxi — 


AttR^ 


roo r 

dz/^a,(W). 


The lower limit of the integral is the threshold ionization frequency of the ion X*. cTj,(X*) is 
the photoelectric absorption differential cross section of the ion X*. 

In this simple scheme the ionization state of the gas is completely determined, at equi¬ 
librium, by the value of the ionization parameter U = Fui/ueC. Q 

With these approximations the ionization parameter U can be written as a function of 
the ratio between any two consecutive ionic species of the generic element X: 


( 5 ) 


U 


\ / Fhi\ (y-rec{X\Te) 
/ \ Fxi ) c 


^Neglecting the normally small effects of Auger ionization, collisional ionization and three body 
recombination. 

"^The inverse of the photoionization and radiative recombination rates are usually referred in the literature 
as the photoionization and recombination times: 


(3), (4) 


fX' _ 
'-1 ~ 


F 


X‘ 


t 


Xi 

2 


1 

arec{Xi,Te)ne' 


Instead in this paper we refer to different quantities as photoionization time and recombination time: see 
eq. 7. 
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2.3.2. Time-evolving photoionization: Equations 


If the gas is not in equilibrium eq. 5 is of course meaningless. However, it is still possible 
to use it formally by introducing an ionization parameter which depends both 

on time and on the specihc ionic species under consideration. 

The time evolution of the relative density of the ion i of the element X, considering only 
radiative recombination and photoionization, is then given by (Krolik and Kriss, 1995): 


( 6 ) 


duxi 

dt 


Fxi + arec{X^ iTf,)ne Uxi + Exi-lUxi-i + Q;rec(^\ • 


The first term on the right of eq. 6, is the destruction rate of the ion X®, both by 
photoionization X® —X®"*"^ and radiative recombination X® —X®“^, while the other two 
terms indicate the formation rate of the ion X® by photoionization of the ion X®“^ and 
radiative recombination of the ion X®"*"^ respectively. 

The solution of eq. 6 is a system of N coupled integral equations in the N unknowns 
nxi, which is analytically solvable only for N=2, with the addition of the charge conservation 
condition. These solutions dehne the time scale teq, that measures the time necessary for the 
gas to reach photoionization equilibrium with the ionizing continuum. This time is given at 
any point of the light curve of the ionizing continuum by the inverse of the destruction rate 
of the ion X* (following Krolik & Kriss 1995). A useful analytical approximation for teq is: 


1 

t OreciXlljIsl'] f Ilx±ti 

V arec{X\Te) / eg ^ 

where eq indicates the equilibrium quantities. During increasing ionization flux phases we 
call photoionization time., tph] during decreasing ionization flux phases we call t^q'^^^^ 

recombination time, tree- These times are generally different from ti and t 2 in eq. 3 and 4. 

Equation 7 shows that the time ^ t dt) necessary for the gas to reach 

equilibrium depends explicitly on the electron density Ue in the cloud, and on the equilibrium 
ratio between two consecutive ionic species calculated at the time t+dt. This is the key result 
of this work, which has major consequences: 


egJ t+dt 


(7) tfg’^'^\t^t + dt) 


^reci,^ ■) ^qTl^ 



11 


1. the time scale on which the gas reaches equilibrium with the ionizing continuum de¬ 
pends on the electron density, even when the continuum increases (Fig. 5); 

2. different ions reach their equilibrium relative abundances at different points of the light 
curve (Fig. 6); 

3. if changes of intensity are not instantaneous (dL/dt< oo), the time behavior of the 
relative abundance of a given ion can be opposite to that of the ionizing source (Fig. 
6 ). 

All of these effetcs have been seen in ionized absorbers in AGN. Using these effects non¬ 
equilibrium photoionization models can strongly constrain the physical state of the absorber. 


2.3.3. Time-evolving Photoionization: Calculations 


We created a program to solve the first order differential equation system of eq. 6 for 
all the ions of the elements H, He, Li, C, N, O. The program uses an iterative method 
(see Gallavotti, 1983), that permits the solution of any system of N first order differential 
equations in the N unknowns of the form x(r) = f(x(r)), (Vr > 0), with the only 
conditions being that fi G and a limited ensemble G exists, such that Xi G G. 

We consider the photoionization from the K-shell of each element, and radiative recombina¬ 
tion to all levels for each ion. 

We use the recombination rates tabulated by Shull and Van Steenberg (1982) for the 
metals. We take the values of the recombination rate of hydrogen from Ferland (1996); we 
get the total recombination rate by summing over levels n=l,20. We calculate the photoion¬ 
ization rate from the K-shell of each ion, carrying out the integrals in eq. 2, using for the 
spectral shape of the ionizing continuum, from the Lyman limit to y-rays, a simple power 
law with a = 1.3 (similar to the observed SED of NGG 4051). We use the photoelectric 
K-shell cross section tabulated by Kallman and Krolik (1995). 

We carry out the calculation of the time-evolving heating-cooling balance, as described 
below. We calculate with GLOUDY (Ferland, 1996) a grid of models for 300 values of U (from 
0.01 to 100), and build the curve U=U(Te), using the technique described in Kallman and 
Krolik (1995). We then interpolate on these, to obtain the initial self-consistent equilibrium 
electron temperature of the gas. The time evolution of the temperature in the cloud is carried 
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out performing an iterative calculation of the time dependent photoionization equations, 
using the dehnition of {t) given above. 

The inputs the program needs are: (a) the initial equilibrium value of the ionization 
parameter U; (b) the spectral shape of the ionizing continuum; (c) the light curve of the 
ionizing continuum; (d) the electron density Ug; (e) the ratio, P, between the source intensity 
and the intensity the source should have in order to produce the observed degree of ionization, 
assuming equilibrium at the beginning of the light curve. 

The output of the program is a list of the relative ionic abundances of the chosen element, 
and the source flux at the illuminated face of the cloud as a function of time. 


2.3.4- Limits of the code 


We neglect photoionization from the L-shell. This affects the computation of the abun¬ 
dance balance when the gas is allowed to recombine to medium-to-low ionization states, with 
the Li-like ions (or lower) being highly populated. However this is not likely to be the case 
for warm absorbers in AGN, for which the mean degree of ionization is usually very high, and 
the most abundant element (C to Ne) are almost equally distributed between He-like and 
fully stripped ions (see Fig. la). Our code can provide accurate results only at medium-high 
ionization regimes, those of interest for the astrophysical problem discussed in this paper. 

We have also neglected Iron ions. This will likely cause an overestimation of the OVII-IX 
and NeIX-XI K-shell photoionization, since all these ions compete for the same photons. By 
comparing the OVII-IX, NeIX-XI equilibrium abundances obtained with our method with 
those obtained with CLOUDY for U in the range 0.1-50 we estimate that the percentage 
differences are smaller than 7 % 

The iterative technique presented in the previous section is only an approximation to 
the correct self-consistent time-evolving heating-cooling calculation accounting for all the 
physical and dynamical heating-cooling mechanisms. We verihed that for high gas electron 
densities (and therefore in a situation close to equilibrium) this method is rather accurate. By 
comparing the equilibrium relative abundances of the main ions obtained with our method 
with that obtained using CLOUDY we estimate that for U in the range 0.1-50 the method 
works with a precision better than 10%. In the opposite case, when the density is sufficiently 
low, the gas remains basically in the same state and does not respond to variation of the 
ionizing continuum. Of course in this case our approximation is also very good, albeit 
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uninteresting. In the intermediate cases, when the heating timescale may be longer than 
the typical photoionization timescale, the recombination rate coefficients may be greater 
than those estimated by our assumption. We have tried to estimate the magnitude of the 
uncertainty on the ion abundance distribution induced by this effect by performing the 
time-evolving calculations hxing the temperature at a constant value equal to the initial 
value, and letting the ionizing intensity be free to vary by a factor of 10. For an average 
gas ionization state typical of AGN warm absorbers (U = 0.1-50), the percentage difference 
between the relative abundances of the main ions calculated with the iterative solution of the 
time dependent photoionization equations, and those calculated for a gas with a constant 
temperature are ^ 15 — 20%. The difference is so small because of the rather small dymamic 
range spanned by the temperature in a photoionized gas with U=0.1-50 (3 x 10^^Te<3 x 10®) 
and of the weak dependence of the recombination rates on Te (about the square root of T^, 
see Shull and Van Steenberg, 1982). 

Given the conhguration of the gas we considered (a thin, plan-parallel slab obscuring 
the line of sight), and the range of electron and column densities, the light-crossing time 
is smaller than both the typical source intensity variability timescales and the associated 
photoionization recombination timescales. We can then neglect the effects of the light- 
crossing time through the cloud of photoionized gas (see Binette, 1988). 


2.3.5. The Step Function Light Curve 


The simplest case is that of a two state light curve. Let us suppose that the ionizing 
continuum intensity goes instantaneously from a “low” to a “high” state and comes back to 
the “low” state after a time t^^ar (Fig. 5a). The time behaviour of the ionization state of the 
gas irradiated by this continuum depends on the value of the ratio {t ^ t + dt), 

and the amplitude of the flux variation (here we adopt a factor 10 change in flux). 

We considered the case of an optically thin cloud of gas with an initial ionization such 
that the most relevant ionic species of the oxygen are OVIII and OIX (corresponding to 
equilibrium value of U^5 with the adopted SED). 

The time behavior of the relative density of noix is shown in Fig. 5b. In both panels 
different lines identify different values of the electron density, Ue = 10^, 5 x 10^, 10®, 10® 
cm“®. In the upper panel different Ue imply four values of the distance of the cloud from the 
ionizing source. 

The photoionization time (eq. 7) of the gas becomes progressively longer as 
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the density decreases: from about 3 x 10^ s for rif. = 10® cm“®, to ~ 10^ s for Ue = 5 x 10^ 
cm“® to even longer timescales for rie = 10^ cm“® (for ne^lO® cm“® the changes of noix 
during the hrt 10^ s are < 10 %). Note that this is not true for ti in eq. 3, which is the 
dehnition of photoionization time usually found in litterature. 

Formally the dependence of tph on Ue is introduced by the choice of a particular set of 
boundary conditions when solving eq. 6, and hence dehning (t ^ t + dt). Physically, 

hxing the boundary conditions of the system means to choose a particular initial ionic dis¬ 
tribution in the gas, and hence its initial ionization state. Different values of Ue, given the 
initial ionization state of the gas (and hence the ratio between the flux of ionizing photons 
at the illuminated face of the cloud and the electron density), imply different distances of 
the gas from the X-ray source. This is clear in the upper panels of Fig. 5 where the ionizing 
flux at the illuminated face of the cloud is plotted for different values of Ue- 

Recombination times are generally longer > tvar), and can be order of mag¬ 

nitude longer. At the highest density tested {rie = 10^ cm“®) the ionization state of the gas 
is able to relax to the initial equilibrium state in less than 10^ s after the source switch off. 
Recombination time scales for Ue < 10® cm“® are long, the order of many times 10^ s. Since 
the source switch off is instantaneous the relative density of OIX never increases after the 
switch off. 

This case illustrates clearly how photoionization recombination timescales can have 
a strong effect on the changes observed in ionized absorbers and why the photoionization 
timescale depends on electron density. In the following section we present a more realistic 
light-curve, and discuss in detail the main features of our models. 


2.3.6. The gradual rise & deeay light curve: dL/dt< oo 


We now consider a more realistic light curve. In this case the source intensity goes from 
a low state to a high state in a hnite time (4,000 s), and after 2,000 s comes back to the 
initial low state with the same absolute gradient (Fig. 6a). The entire up & down cycle lasts 
10,000 s. As in the previous case the change in flux is a factor 10. The corresponding light 
curves of the relative abundances of the fully stripped ions of three different elements, CVII, 
NVIII and OIX, are shown in Fig. 6b. In both panels different lines correspond to different 
values of the electron density, = 10®, 10® cm“®. 

The gradual changes of the ionizing continuum produce time delays between the source 
light curve and the relative ion abundance light curves. In the lower density case (ue = 10® 
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cm“^ solid lines), uqix reaches its maximum value (well below the equilibrium value ~ 1) 
around the minimum of the luminosity intensity (at ~ 10, 000 s). During the whole decreasing 
luminosity phase npix is slightly increasing or constant. X-ray spectra taken during the high 
luminosity phase and the decreasing luminosity phase would show an OVIII absorption edge 
correlated with the intensity of the ionizing continuum, as seen in MGC-6-30-15 (Fabian 
et al., 1994, Reynolds et al., 1995). At this density CVII is able to reach its maximum 
equilibrium value (Fig. 6b), but even this ion does not relax to its initial equilibrium value 
for many times up and down cycle time. Fig. 6b also shows that the ions reach their 
maximum values at different times, ricvii reach its maximum value about 4,000 s before 
noix- 


At higher densities (rze = 10® cm“^, dashed lines), the fractional abundance of each of 
the three ions reaches its maximum equilibrium value during the hrst 6,000 s, and relaxes 
to its minimum equilibrium value during the following 2.4 x 10^ s, but with different times- 
scales: OIX reaches its minimum equilibrium value after a time corresponding to two cycles, 
about one cycle after CVII. This could help to explain why models in which all ions of all 
elements are in photoionization equilibrium, so often fail to describe AGN spectral evolution. 
Spectra accumulated immediately after a very steep decreasing intensity phase could contain 
no signihcant absorption features at the energies of OVII-OVIII K-edge (the oxygen being 
completely ionized), but still show a deep absorption edge at E~ 0.5 keV, due to the presence 
of a large amount of recombined CVI-CVII in the absorbing gas. High quality spectra would 
allow powerful tests of non-equilibrium photoionization models. 

A general result is that the observation of any delays in the response of the absorber 
to flux changes on time scales of ~ 5000 — 1000 s, immediately implies photoionization with 
a density in a reasonably restricted range, 10® tUet 10® cm“^ (depending on the average 
ionization state). 


3. Modeling the ROSAT data of NGC4051 


As an example of applicability of our models we present here the case of a ROSAT 
PSPC observation of the low luminosity, rapidly variable Seyfert 1 galaxy NGC4051. 

Both the soft and hard X-ray flux of NGC4051 vary by large factors (up to 20) on a time 
scale of hours (e.g. Lawrence et ah, 1985, Guainazzi et al., 1996) in a roughly correlated way. 
The presence of an ionized absorber in this source was hrst proposed on the basis of variations 
of the GINGA softness ratio correlated with the hux (Fiore et al. 1992). ROSAT PSPG 
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observations, and snbseqnently observations with the higher resolntion CCDs on board of 
the ASCA satellite (Ponnds et ah 1994, Me Hardy et ah 1995, Mihara et al. 1994, and 
Gnainazzi et al. 1996), also snggested the presence of a high colnmn density, highly ionized 
absorber throngh the detection of a deep absorption edge at 0.8-0.9 keV. 

McHardy et al. (1995), hnd that a simple single zone eqnilibrinm photoionization model 
can provide a reasonably good representation of spectra accumulated in 1000-3000 seconds, 
but the best £t ionization parameter does not track the source intensity, as required by the 
model. This behaviour cannot be explained in a single zone ionization equilibrium model and 
therefore NGC4051 is a good target to test our time-evolving photoionization and collisional 
ionization models. We decided to compare our models to the ROSAT PSPC data aquired 
on 1991 November 16, and reported by McHardy et al. (1995), when the source count rate 
showed large and rapid variations (up to a factor of 6 in a few thousand sec., see figure 5 of 
McHardy et al., 1995). Analysis of the ASCA observation is deferred to a future pubblication. 

The data reduction and the timing analysis were performed using the PROS package in 
IRAF. The observation spanned 77 ksec and contained 28.7 ksec of exposure time. NGC 4051 
gave a mean count rate of 1.6 s“^. We accumulated eight spectra {a — h), using a 3’ radius 
extraction region, accumulating contiguous data with similar count rates. Background counts 
and spectra were accumulated from an annulus of internal and external radius of 3’.5 and 6’ 
respectively. 


3.1. Hardness ratio analysis: the absorber is not in photoionization 

equilibrinm 


Independent of any spectral fit the behavior of the main physical properties of the 
absorber can be seen in a color-color diagram. In Fig. 7 we plot the hardness ratios HR=H/M 
against the softness ratio SR=S/M from the count rates in three bands (S=0.1-0.6 keV, 
M=0.9-1.5 keV, and H=1.7-2.5 keV) for theoretical curves (for log(Ah/)=22, 22.5 and 23) 
obtained by folding the equilibrium photoionization models (for log(U) in the range —0.3 
to -1-1.5, and Galactic Nh'- 1.31x10^° cm“^, Elvis et al. 1989) with the response matrix of 
the PSPG. We also plot the colors of the source in the eight spectra (a - h). All the data 
points he in a region of this diagram corresponding to the high U ends of the photoionization 
theoretical curves where both SR and HR decrease, as U increases, until all the ions in the gas 
are full stripped and the gas is completely transparent to radiation of any energy. The colors 
of NGG 4051 are all consistent with the OVHI-OIX ions being dominant. Filled circles on 
the logA’//=22.5 curve mark values of U in the range 4.0-7. All the data points are between 
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the two extreme values of U, a factor ~ 1.5 change, while the source intensity varies up to a 
factor ~ 6 in the eight spectra. The gas is clearly not responding to the continuum variations, 
a conclusion equivalent to that obtained by McHardy et al. (1995) using spectral hts with 
equilibrium photoionization models (see their hgure 5). In the framework of the models 
discussed in this paper this behavior suggests three different possibilities (we do not take 
in account pure collisional ionization in this case because the observed spectral variations 
should be attributed to ‘ad-hoc’ changes of Te or on time scales as short as 2000-4000 
s): 

1. the gas is far from ionization equilibrium; 

2. the gas has a distribution of densities; 

3. both collisional and photoionization processes are comparably important in the same 
physical region. 

We investigated these possibilities in turn using detailed spectral fits. 


3.2. Time-evolving photoionization 


From Fig. 1 we see that a given ionization state can be roughly determined by the 
measure of at least two consecutive ion abundances, e.g. OVII and OVIII. The measure of a 
single edge in fact would not distinguish between ‘low’ and ‘high’ ionization solutions. The 
same measured feature could be produced by a lower Nh, lower mean ionization gas, or by 
an higher Nh, higher mean ionization state gas. 

The best derived quantities to compare observed spectra with our several physical mod¬ 
els are atomic edge strengths. Here we are mainly interested in the OVII and OVIII because 
they are the strongest and therefore the easiest to detect and measure. However, with the 
PSPC OVII and OVIII edges are not individually discernable, and we must resort to model 
hts with multiple components. Onr choice is to use the components that can ensure an 
estimate of the OVII and OVIII r as nnbiased as possible. 


3.2.1. A 3-edge ‘Model Independent’ Speetral Fit 
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We fitted the eight spectra with a model consisting of a power law rednced at low energy 
by the Galactic colnmn density, the OVII and OVIII edges and another edge at 1.36 keV 
to acconnt for possible spectral complexity in the 1-2 keV region (in particular the NeIX-X 
K and Fe XV-XX Lp absorption discussed in §2.1). Five parameters were at first allowed 
to be free to vary: the spectral index the OVII and OVIII edges r, the 1.36 keV edge 
r and the model normalization. The results are presented in Table 3. The hts with the 3 
edges model are acceptable in all cases. We stress that r(1.36 keV) in Table 3 should not 
be regarded as a true measure of the optical depth of Ne K and Fe L ions. This feature 
provides only one of the possible parameterizations of the spectrum in the 1-2 keV range, a 
band in which a change in the continuum spectral index could also be present (the 2-10 keV 
Ginga and ASGA spectra of this source are typically flatter by ~ 0.5 than the PSPG 0.1-2 
keV spectra). As discussed in §2.1 is difficult to discriminate between Ne absorption and a 
real spectral flattening above 1 keV with instruments of moderate spectral resolution like 
the PSPG. We have performed a series of fit using a broken power law with break energy in 
the 1-2 keV band and the two oxygen edges, obtaining r similar to those reported in Table 
3. We are therefore conhdent that the estimation of the oxygen edge r is robust, within the 
rather large uncertainties given in Table 3. In principle, a way to reduce the uncertainties 
is to £x the continuum spectral index to a common value. The results of this series of fits 
are again given in Table 3. The uncertainies on r(OVII) and r(1.36 keV) are indeed smaller 
than in the previous case but this is not true for the r(OVIII) uncertainties. The reason is 
that there is a strong anti-correlation between aE and tovii- A similar anti-correlation is 
present between aE and r(1.36/ceV). In contrast no correlation is present between aE and 
r(OVIII). This is illustrated in Fig. 8 where we show the contours of these parameters 
for spectrum g. Since the r(OVII) values could be biased in the fit with hxed q;_b by the 
r(OVII)-Q;£; correlation, we prefer to compare our time-evolving models to the OVII and 
OVIII optical depths obtained leaving aE free to vary. 

In all but one case (spectrum d) the depth of the OVIII edge is higher or comparable to 
that of OVII, suggesting a ‘high’ ionization solution (U>4, see Fig. 1), consistent with the 
hardness ratio analysis of Fig. 7. 


3.2.2. Comparison between models and the oxygen edge depths: evidenee for a 

non-eguilibrium photoionization absorber? 


We converted the best fit r into OVII and OVIII relative abundances assuming a solar 
oxygen abundance and a total hydrogen column density Nh- An indication of Nh comes from 
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Table 1; NGC 4051: 3-Edge Model Fits 


Spectrum 

r(0.74 keV) 

r(0.87 keV) 

r(1.36 keV) 

CtE 

X^(d.o.f.) 

a 

U.D_o.4 

A 04+0.29 
^•’^^-0.23 

O.S+JJ;^ 

0.84l°;|i 

1.0±0.4 

0.81++ 

1.3±0.1 
1.34 (fixed) 

0.84(21) 

0.82(22) 

b 

< 0.5 

0.66±0.21 

0.66±0.21 

0.7±0.3 

0 71+01® 

‘ 1-0.16 

1.3±0.1 
1.34 (fixed) 

1.19(23) 

1.14(24) 

c 

U.O_o.4 

U.00-0.26 

A +0.31 

U- + -L_0.34 

A 4-1 +0.31 

U.^1-0.34 

< 0.6 

< 0.42 

1.3±0.1 
1.34 (fixed) 

0.85(22) 

0.82(23) 

d 

< 0.7 
o.99t°J2 

< 0.3 

< 0.32 

1 1+0.8 
1-1-0.6 

1 Q^+0-®l 
i.yo_0.54 

1-1-0.2 
1.34 (fixed) 

0.95(18) 

1.0809) 

e 

U.OZ_o 22 

0.86±0.28 

0-86l°j® 

0.5±0.3 

0.65±0.19 

1.4±0.1 
1.34 (fixed) 

1.30(22) 

1.2503) 

f 

< 0.6 

0 f:4+0'20 

1-01+J^ 

1.UZ_o.25 

< 0.5 

0.44±0.14 

1.4±0.1 
1.34 (fixed) 

1.41(23) 

1.42(24) 

g 

< 0.5 

< 0.18 

1.47±0.29 

1-42+;24 

< 0.5 

< 0.21 

1.3±0.1 
1.34 (fixed) 

0.65(22) 

0.6703) 

h 

1 a+0.6 

-‘-•^-0.5 

0.63t°;i 

0.79l°l 

0-86l°;i 

0.9+;^ 

0.45+:?o 

1.2±0.1 
1.34 (fixed) 

0.72(22) 

0.8203) 
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the color-color diagram of Fig. 7. Although calculated using a photoionization equilibrium 
model, the theoretical curves in this diagram suggest a value for logN^ between 22 and 23, 
and so we adopt logA^i^=22.5. The three panels of Fig. 9 show the light curves of the source 
count rate (upper panel) and of the OVIII and OVII abundances (middle and lower panels 
respectively). 

The time evolution of the ionization structure of a cloud of gas photoionized by a variable 
source is complex and its behaviour sometime counter-intuitive. We then discuss hrst the 
simplest case: high electron density, for which each ion is close to its equilibrium state. We 
examine next the case of lower densities and hence non-eqnilibrium solutions. 

The dotted curves on the middle and lower panels of Fig. 9 show the Ug = 10^'^ cm“^, 
P=1 (the ratio of the incident flux to that needed to produce the initial ionization distribntion 
assnming photoionization eqnilibrinm) curves, when the gas is close to eqnilibrinm. While 
the equilibrium OVII curve tracks the count rate variations (it is strictly anticorrelated with 
the count rate light curve), the OVIII curve does not. The different OVIII behavionr is due 
to the different balance in the destrnction rates of OIX and OVIII. When the ionizing flux 
is at its maximum most of the oxygen is OIX. When the flnx decreases from the maximnm 
(from point b to d) OIX recombines to OVIII increasing the OVIII abnndance. When the 
flux decreases from point e to point / at first OVIII increases again because of the high 
destruction rate of OIX, but after a certain point the amonnt of OVIII recombining to OVII 
start to be higher than the amonnt of OIX recombining to OVIII, and so the total OVIII 
abundance start to decrease. Instead, the amonnt of OVIII recombining to OVII is always 
higher than the amonnt of OVII recombining to OVI. This gives rise to a different behavionr 
of the OVIII and OVII curves in response to the same ionizing flux variations. It is interesting 
to note that in this case while the dynamical range of variation of the OVII curve is larger 
than that of the ionizing flux, the OVIII equilibrium curve shows a more compressed range 
of variations (when the ionizing flux varies by a factor of 6 the OVIII and OVII abundances 
vary by a factor of 3 and 30 respectively, see Fig. 9). Therefore variations of OVIII would be 
much more difficult to detect than variations of OVII, at least at these regimes of ionization. 
We re-emphasize: the behaviour of a single edge does not provide a unigue interpretation of 
the data. 

As explained in §2.3.3, a grid of theoretical OVII and OVIII light curves was generated 
using our time-evolving photoionization code for 28 valnes of Ue from 5 x 10® cm“^ to 10® 
cm“^ and 15 valnes of P from 0.5 to 2. 

We compared these curves with the measured relative abundances of OVII and OVIII and 
found the best fit model using a technique. The thick solid lines in the lower and middle 
panels of Fig. 9 represent the best fit non-equilibrium OVII and OVIII cnrves. The agreement 
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between the best fit model and the observed OVII and OVIII abundances is good: = 1.17 

for 14 dof. The best htting values for Ue and P are tightly constrained: Ue = (I.OIqJ) x 10® 
cm“®, P= 1.5lo;3. From the best £t rie we can estimate the distance of the absorbing cloud 
from the central source. We obtain R= (0.74i'io;4o) x 10^® cm (3 light days). 

The dashed lines represent the solutions obtained using the 1 a confidence interval on 
TT-e and P. The best fit curves (and the 1 a confidence intervals curves) show a compressed 
dynamical range of OVIII and OVII abundances variations and a delay between the source 
maxima and the ion abundance minima of 3000-6000 sec. The compressed dynamical range 
is due to a mean over-ionization of the gas. While the best fit P shows that the gas in the 
initial point a is near to equilibrium, it departs strongly from equilibrium during the later 
low intensity states (spectra d to h). So, despite the fact that the source spends more time 
in low states than in high states, the gas density is sufficiently low that the gas does not have 
time to fully recombine after the few events when it suffers high illumination and becomes 
highly ionized. 

The above results were obtained assuming a total hydrogen column density of log(iVji^)=22.5. 
Assuming an higher (lower) column would imply a mean lower (higher) OVII and OVIII rel¬ 
ative ion abundance. Therefore, in principle the accurate measure of both edges would 
constrain also the total warm column density. The uncertainties on the PSPC determi¬ 
nations however preclude this possibility, and better resolution measurements are therefore 
needed. The energy resolution of the ASCA SIS, for example, is just sufficient to separate the 
OVII and OVIII absorption edges. A quantitative test of non-equilibrium photoionization 
model using the ASCA data and a comparison between ROSAT PSPC and ASCA data of 
NGC 4051 is beyond the scope of this paper and will be presented in a future pubblication. 

A much better separation, and therefore characterization, of the absorption features will be 
possible with the high resolution (factor of 10-30 better than ASCA SIS) gratings on AXAF 
and XMM. 


3.3. Other models 


Despite the success of time-evolving photoionization models, alternatives do exist. We 
discuss two of them in the following. 


3.3.1. Large density variations in the absorbing gas 
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If the gas is not confined to a single cloud with constant density but rather is distributed 
in a region with, say, an increasing density, then different ionization equilibria could apply 
to different regions in the cloud. Two extreme regions may exist: in the region with higher 
density, lower U, collisional ionization will be the dominant ionization mechanism and the 
spectra transmitted by this region would show always the same features, irrespective of the 
source intensity; the other region, with lower density and higher U, is completely ionized (for 
carbon, oxygen, neon and iron ions up to FeXXII) when the source is in the high state, but 
when the source is in a low state the abundances of OVIII are sizeble and imprint the edges 
in the spectrum seen in low intensity spectra. Here we are assuming that the density of the 
photoionized region is high enough for the gas to be in instantaneous equilibrium with the 
ionizing intensity. If not, the average ionization degree of the gas would be very high during 
the whole observation, and the gas would be always almost transparent at the energies of 
the relevant absorption edges. 

We tested this hypothesis by htting the highest intensity spectrum (6) with a simple 
power law model plus a collisionally ionized absorber, using the method of Fiore et al. (1993). 
The best £t temperature and are 2.8 ± 0.1 x 10® K and l.Slg;^ x lO^^cm”^ respectively. 
We then used the same model (with fixed continuum parameters, fixed temperature and 
Nh but variable normalization) with the inclusion of an additional OVIII edge to mimic a 
variable ionization state, equilibrium photoionization absorber. The fits are all acceptable. 
The worst fit is that of spectrum (d) (x2 = 1.4, 21 dof, probability of 10.4 %). In Fig. 10 
we plot t{OVIII) as a function of the time. The dynamic range of variations of r is here 
larger than that on noviii in Fig. 9 but is still more compressed than that predicted by 
equilibrium photoionization model (solid line). We can therefore exclude the possibility that 
a major part of this absorber could be in pure photoionization equilibrium with the ionizing 
intensity. 


3.4. A “hot” photoionized absorber 


The other possibility is that both collisional and photoionization processes are important 
in the same physical region. In this case the transmitted spectra are very similar to those 
transmitted from purely photoionized clouds of gas which much higher ionization parameter 
(§2.1, Fig. 4). The electron temperature of the gas is mainly determined by collisions, and 
is higher than that expected in pure photoionization equilibrium. The ionization parameter 
U is no longer linearly correlated with the ionizing intensity, and then its dynamical range 
of variations is compressed by a factor > 2 compared to the pure photoionization case. 
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This could in principle explain the compression observed in the measured dynamical 
range of variation of noviiii but could not account for the delays observed on the response 
of the ionization state of the gas to source intensity variation. Unfortunately there is no way 
to distinguish between ‘hot’ photoionized absorber models and simple pure photoionization 
models on the only basis of the spectral analysis. Nevertheless the delay observed in the 
response of the ionization degree of the gas between spectrum (b) and (c) suggests that a 
non-equilibrium photoionization component is mainly required by data. 


4. Conclusion 


We have investigated ionization models for AGN in different regimes of gas volume 
densities and photoionization states. In particular we focussed on ’high gas density, low pho¬ 
toionization parameter’ gas clouds, where collisional ionization is likely to play a signihcant 
role in the gas ionization, and on low gas densities, where the photoionization may be far 
from equilibrium. 

We presented detailed model calculations in both regimes. While the time-evolving pho¬ 
toionization models in §2.3 are far from being complete or exhaustive they are nevertheless 
instructive, and reveal the main features of these kinds of models. 

Our main hndings can be summarized as follows: 

1. In many important astrophysical conditions (OVII, OVIII regime) the fractional abun¬ 
dances of the most important ions of O and Ne in photoionized gas are more broadly 
distributed in U than those of collisionally ionized gas are in T. 

2. In the collisional ionization case the heavy elements show a strong inertia against 
becoming highly ionized, even when lighter elements, like Oxygen, are almost fully 
stripped. In this case the transmitted spectrum shows a large and complicated absorp¬ 
tion structure between 1 and 3 keV, mainly due to iron L XVII-XIX and Neon K IX-X 
absorption, which is much less visible in spectra emerging from photoionized gas with 
similar OVII and OVIII abundances. This absorption structure may be mis-interpreted 
as a flattening of the spectra above ~ 1 keV, when fitting low energy resolution data 
with a photoionization equilibrium model. Higher spectral resolution and good S/N 
observations are therefore needed to distinguish between collisional ionization and pho¬ 
toionization. 
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3. In non-equilibrium photoionization models the response of the ionization state of the 
gas to sudden changes of the ionizing continuum is delayed even during increasing 
luminosity phases. The delays increase for decreasing electron densities, as changes of 
He, require changes of the intensity of the ionizing flux, i.e. changes of the distance of 
the gas from the X-ray source (taking as hxed the initial ionization state of the gas). 

4. The recombination timescale is generally much longer (up to orders of magnitudes) 

than the photoionization timescale, because of the dependence of {t ^ t + dt) 

on the equilibrium ratio {nxi+^/nxi), evaluated at the time t + dt. This means that a 
photoionized absorber undergoing frequent, quick, and consistent changes of ionizing 
luminosity is likely to be overionized with respect to the equilibrium ionization, a state 
that would be reached only after a sufficiently long low intensity phase. 

5. If the changes of the ionizing luminosity are not instantaneous, and the electron density 
of the cloud is low enough, the ionization state of the gas could continue to increase 
during decreasing source luminosity phases. This means that we may measure a max¬ 
imum in ionization state of a given element, when the ionizing flux is at a minimum 
(opposite to what is expected in equilibrium models). 

6. Different ions of different elements reach their equilibrium abundance on different 

timescales. This is again because of the dependence of {t t + dt) on the 

ratio {rixi+i/nxi)■ Therefore in the same cloud of gas carbon could be in equilibrium 
while oxygen could be very far from equilibrium. This may help in explaining why 
models where all ions of all elements are in photoionization equilibrium so often fails 
to provide a reasonable description of AGN spectra and spectra evolution 

We have tested the above models in the case of the Seyfert 1 galaxy NGC4051. The 
ROSAT observations of NGG4051 are not consistent with a simple equilibrium model, but 
can be explained straightforwardly by our time-evolving photoionization models. The two 
main features in the non-equilibrium best £t models are: (a) the compressed range of vari¬ 
ability of the measured OVII and OVIII relative abundances with respect to the amplitude 
of the source variations, and to the amplitude of the variations of the abundances of these 
ions expected in equilibrium photoionization models; (b) the 3000-6000 sec delay between 
the maximum intensity state of the source (spectrum b) and the minimum of the best fit 
OVIII abundance curve, i.e. the maximum ionization state of the gas (spectrum c). As 
result we were able to estimate the gas electron density, = (l.OtgJ) x 10® cm“® (assum¬ 
ing logNH=22.5) and hence the distance of the ionized gas cloud from the X-ray source in 
R= (0.74tQ;4Q) X 10^® cm (3 light days). We explored alternative models and we also ex¬ 
plored alternatives which we find to be less likely; we discuss ways to distinguish between 
them conclusively. 
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We conclude that non-equilibrium photoionization and collisional models apply to wide 
zones of gas density and ionization, zones which are expected in AGN. These effects must be 
considered in understanding ionized absorbers, and seem likely to explain otherwise puzzling 
behavior, without resorting to ad hoc distributions of gas. Several clear diagnostics of these 
models exist so that decisive tests will soon be possible. 
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Fig. 1.— Fractional abundances of NeVIII-NeXI, and OVFOVIX, calculated in the case of 
photoionization (Fig. la: upper panel), and collisional ionization (Fig. lb: lower panel). 
The two intervals of U and Te highlighted are: (a) the interval for which both OVII and 
OVIII abundances are greater than 0.2 (dotted lines), and (b) the one for which the OIX 
relative abundance is greater than 0.75. 
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Fig. 3.— (a) The ratio between two power law spectra emerging from clouds of collisional 
ionized (Tg = 1.8 x 10® K) and photoionized (U=2) gas with solar abundances, (b) The 
ratio between two power law spectra emerging from clouds of collisionally ionized gas with 
Tg = 3.9 X 10® K, and photoionized gas with U=10. The column density is Nh = 10^^ ® cm“^ 
in both cases. 
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Fig. 4.— (a) Spectra from gas in which both collisional ionization and photoionization are 
important. The temperature of the gas is hxed at = 2.5 x 10® K, and log(U)=-l, -1, 0.3, 
0. 0.4. (b) Spectra from gas in pure photoionization equilibrium., with log(U)=0.7, 1.1, 1.4. 
The column density is Nh = 10^^ ® cm“^ in both cases. 
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Fig. 5.— Upper panel; step light curve of the ionizing radiation, for four different values of 
the electron density (ng = 10^, 5 x 10^ 10®, and 10® cni“®, which correspond to four different 
values of the incident flux). Lower panel: time behaviour of the relative abundance of OIX, 
for the four values of rig. 
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Fig. 6.— Upper panel: 2-phases light curve of the ionizing continuum. Lower panel: time 
behaviour of the relative abundances of the fully stripped ions of three elements: CVII, 
NVIII, OIX. In both panel different lines correspond to different values of the electron density, 
rie = 10® (solid line), and 10® cm“® (dashed line). 
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Fig. 7.— Color-color diagram of the eight spectra of NGC4051. Lines are theoretical pho¬ 
toionization curves, built folding equilibrium photoionization models with the response ma¬ 
trix of the PSPC. Different lines correspond to three different values of the gas column den¬ 
sity: log(7V//)=22.5 (solid line), log(iV//)=22 (dotted line), and log(iVH)=23 (dashed line). 
On each curve the ionization parameter U, increases going from top to bottom. U values are 
indicated on the log(iVj:^)=22.5 line. 
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Fig. 9.— Light curve of NGC 4051 during the November 1991 ROSAT-PSPC observation 
(Upper panel). Middle and lower panel; Best £t noviii and novii abundances (points 
with errorbars), obtained by htting a power law continuum plus 3-edge model to the 8 
spectra, and by using log(A'H)=22.5. Solid lines in the two panels represent the best £t non¬ 
equilibrium photoionization model to the noviii and novii data respectively, while dashed 
lines are the solutions obtained using the 1 a conhdence intervale on Ue e P. Dotted lines 
are the equilibrium noviii and novii curves, and was obtained using the non-equilibrium 
photoionization code by hxing rig = 10^° cm“^ and P=l. 
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Fig. 10.— Similar to Fig. 9 but the lower panel shows the OVIII edge r obtained htting 
the eight spectra with a model consisting of a power law with energy index 1.3 (fixed in the 
fit) plus a collisional ionized absorber with Nh = 1-5 x 10^^ cm“^ and T=2.8 x 10® K (fixed 
in the £t) plus an edge at 0.87 keV. Two parameters only are allowed to vary: the model 
normalization and the OVIII edge r. 




















